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Abstract. Abundance analysis of carbon has been performed 
in a sample of 80 late F and early G type dwarf stars in 
the me tallicity range —1.06 < [Fe/H] < 0.26 using the forbid- 
den [C I] line at 8727A. This line is presumably less sensi- 
tive to temperature, atmospheric structure and departures from 
LTE than alternative carbon criteria. We find that [C/Fe] de- 
creases slowly with increasing [Fe/H] with an overall slope 
of —0.17 ± 0.03. Our results are consistent with carbon en- 
richment by superwinds of metal-rich massive stars but incon- 
sistent with a main origin of carbon in low-mass stars. This 
follows in particular from a comparison between the relation 
of [C/O] with metallicity for the Galactic stars and the cor- 
responding relation observed for dwarf irregular galaxies. The 
significance of intermediate-mass stars for the production of 
carbon in the Galaxy is still somewhat unclear. 
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1. Introduction 

The stellar origin of carbon became clear with the discovery of 
the Triple Alpha reaction (Opik 1951; Salpeter 1952) and the 



beautiful demonstration by Hoyle (1954) that the high carbon 



abundance would require a resonance state at 7.65MeV of car- 
bon, a prediction that was soon verified experimentally. But in 
which type of stars was carbon formed? 



In their classical paper Burbidge et al. ( 1957 ) suggested 
carbon to be provided by ma ss los s from red giants and su- 
pergiants. Arnett & Schramm ( 1973 ) argued that massive stars 
might be most effective. They discussed the composition of 
matter estimated to be ejected from stars (with helium core 
masses above 4A4©, i.e. A^tot > WMq) and found C/O ra- 
tios consistent with the solar system values. The yield prescrip- 
tions developed by Talbot & Arnett (1973, 1974) consequently 
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ascribed carbon production to the massive stars. However, Tru- 



ran ( 1 977 ) and Tinsley ( 1 977 ) argued that carbon stars of lower 



mass could also be significant sites for carbon. Dearborn et al. 



(1978) also suggested that low mass stars (as evidenced by car- 
bon stars and carbon-rich planetary nebulae) may be a signif- 



icant source of 12 C, while Iben & Truran ( 1978 ) concluded 
from thermally pulsing models that intermediate mass stars and 
high mass stars contributed carbon in roughly equal amounts. 
The arguments for significant contributions of carbon by low- 
mass an d inte rmediate-mass stars were fur ther e laborated by 
Tin sley (1978 ) and Sarmiento & Peimbert ( 1985 ). Timmes et 
al. ( 1995 ) in their models of Galactic chemical evolution with 
then available yields for supernovae as well as for low- and 
intermediate mass stars, found very significant contributions 
of carbon from the latter to the Galactic Disk. (For a general 
review on the significance of low- and intermediate-mass car- 
bon stars for Galactic nucleosynthesis, see G ustafss on & Ryde 
(1998).) Recently, Kobulnicky & Skillman fll998| ) concluded 
from the correlation between C/O and N/O ratios found in H 
II regions in three starforming galaxies, that the formation of 
N and C is coupled and presumably occurs in intermediate or 
low mass stars. Prantzos et al. (1994), however, adopting the 
idea of Maeder (1992) that radiatively driven massive winds 
from high-mass stars should provide huge amounts of helium 
and carbon, found these stars to be the main contributors while 
the role of intermediate-mass stars seemed much less signifi- 
cant. Similarly, Garnett et al. (1995) observed that the C/O ra- 
tio increased with increasing O/H for dwarf irregular galaxies, 
a trend which they found to be consistent with the suggestion 
that carbon is produced in massive stars with yields dependent 
on metallicity, as is expected for stars with radiatively driven 
winds. 

The shifting views as regards the relative role of high-mass, 
intermediate-mass and low-mass stars have reflected the differ- 
ent uncertainties and conjectures concerning the carbon pro- 
duction, dredge-up and mass loss from stars of different mass 
and metallicity. Although considerable progress has been made 
in these respects, there are still very severe remaining uncer- 
tainties. Also as regards observed carbon abundances for stars 
and planetary nebulae the uncertainties have been, and are still, 
considerable which makes checks of the theoretical results, as 
well as more direct empirical approaches, difficult. We shall 
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here present an attempt to improve in this latter respect, as re- 
gards carbon abundances for solar-type stars in the Galactic 
Disk. The idea is to use high S/N observations of a hitherto 
little used forbidden carbon line in the infrared as the single 
criterion. The strength of this line, in comparison with other 
carbon abundance criteria, is not at all as sensitive to other pa- 
rameters as to the carbon abundance. We shall also see that 
with high quality data it is possible to carry out more direct 
tests concerning the role of stars of different masses; tests that 
are independent on the results of the calculations of yields. 

2. Observations and data reductions 

The observations were made at European Southern Observa- 
tory (ESO) with the 1.4m Coude Auxiliary Telescope (CAT) 
coupled to the Coude Echelle Spectrometer (CES) and the 
Long Camera during 3 observing runs with the following CCD 
detectors: July 2-15 1994 CCD#30, Dec 29 1994-January 4 
1995 CCD#34 and Sept 14-26 1995 CCD#34. The slit-viewing 
auto-guiding system ensured efficient light collection. The fo- 
cusing and alignment of the CCD cameras were checked before 
each observing night. 

Spectra centered near the 8727A [C I] line and with a 
wavelength range of about 70A and a resolution of 65,000 
were obtained for 87 southern stars of the Edvardsson et al. 



1.05 



(|1993j) (cited as EAGLNT in the following) survey with dif- 
ferent metallicities. Depending on the stellar magnitudes and 
weather conditions integrations of between 20 seconds (Pro- 
cyon, V = 0.37) and 8 times 60 minutes (HD 199289, V = 
8.29) were obtained for the stars. The long integration times 
were motivated by the need to obtain high S/N ratios (prefer- 
ably above 200) for the generally very weak [C I] line. In order 
to increase the S/N ratio, slightly different central wavelengths 
were chosen for separate observations of the same star. A typ- 
ical S/N ratio for the final, summed up spectra was generally 
between 200 and 300. No spectrum had a S/N ratio below 150. 

The data reductions were performed either on line with the 
ESO IHAP software or later with the ESO MIDAS package, 
and contained the following steps: 

- Subtraction of bias and dark signal from the stellar frames 
and of bias from the flat-field frames. 

- Division of the stellar frame by a flat-field frame obtained 
with unaltered set-up, and within hours of the stellar frame. 

- Summation of the relevant CCD rows to produce a one- 
dimensional spectrum. 

- Wavelength calibration by means of a Th/Ar lamp spectrum 
obtained in connection with the stellar frame. 

- A 3rd order spline function was fitted to the pre-defined 
continuum points in the spectrum, and finally the spectrum 
was divided by this function to yield a relative flux scale 
with continuum level of 1.00. 

For stars where more than one spectrum of the region was 
obtained, both the individual spectra and also the rebinned and 
summed total spectrum were alternatively used in the data anal- 
ysis. 
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Fig. 1. The observed and synthetic spectra of HR 1294 (HD 
26491) with the carbon, iron and silicon lines indicated 



3. Analysis 

3.1. Model atmospheres 

The stellar spectra were synthesised with the Uppsala spectrum 
synthesis program, using the opacity-sampling MARCS model 
atmospheres described in EAGLNT. These are line-blanketed, 
one-dimensional atmospheres with mixing-length convection, 
produced under the assumption of LTE. An example of a com- 
parison between observed and synthetic spectra is shown in 
Fig. 1. 

3.2. Atmospheric parameters 

Originally, the atmospheric parameters T e s, \ogg and [Fc/H] 
of EAGLNT were adopted, and £ t was calculated with the for- 
mula given in EAGLNT 

& = 1.25 + 8-10" 4 (T off -6000)-1.3(log 5 -4.5)kms- 1 .(l) 

In EAGLNT, the surface gravities (g p hot) were based on the 
ci index of the Stromgren photometry. However, surface grav- 
ities were alternatively estimated for the stars from the paral- 
laxes given in the HIPPARCOS catalogue ( |ESA, 1997| ) using 
the equation 



logger = 4 log pf + log 



M 

Mr. 



2 log 7T 



off ""0 

0A{V + BC - 0.26) + 4.44. 



(2) 



The masses for the stars were estimated from evolutionary 



tracks by Schaller et al. (1992) and Claret & Gimenez (1992) 



Due to the uncertainty in the models and the few tabulated 
steps in Z, the uncertainty in the mass estimate may be AM = 
O.2A^ . Ng & Bertelli (|1998[) have also estimated masses for 



these stars showing a good agreement compared to our val- 
ues. The bolometric corrections were estimated according to 



Allen (1973 ). For our stars the bolometric corrections are small. 



The V-magnitudes from the HIPPARCOS catalogue have been 
used. 
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Fig. 2. Comparison between log<? p hot given in Edvardsson et 
al. (1993) (EAGLNT) and logg,,. estimated from Hipparcos 
parallaxes. The plusses represent stars with [Fc/H] < —0.6. 
The dashed line indicates a one-to-one relation to guide the eye 

A comparison between and g p hot is shown in Fig. ||. In 
Tabled the differences between logg p hot and log g v are given. 
It is seen that <? p hot is greater in most cases, the mean difference 
is 0.06 dex, with a formal scatter of 0.12 dex. In Fig. |^ one can 
also see that all stars with [Fe/H] < —0.6 have a lower g v than 
Sphot- The systematic trend that low logg values become even 
lower when estimated from parallaxes might be due to errors in 
the calibration of Stromgren photometry. 

We have here chosen to subsequently use logg p hot since 
the statistical errors in logg^ are greater. However, systematic 
errors in the surface gravities may affect the carbon abundance 
determinations and the comparison with log g^ could give a 
hint as regards the effects of systematic errors. This is discussed 
further below. 

3.3. Atomic and molecular line data 

The [C I] 8727A line is a weak feature, located in the far wing 
of a stronger Si I line, and care must be exercized so that all 
relevant blends affecting the line as well as the neighbour- 
ing continuum regions are included in the synthetic spectra 
used in the analysis. The line list used contains all lines ex- 
isting in the VALD data base (Piskunov et al. 1995) between 
8705A and 8750A. In addition, data for more than one hundred 
CN lines, in the interval 8722A to 8732A, were added to the 
list by Bertrand Plez. Theoretically calculated wavelengths of 
molecular lines are difficult to obtain to high accuracy. Labo- 
ratory measured wavelengths were found for three of the CN 
lines in this interval while the rest were corrected in a semi- 
empirical way using a program that Patrick de Laverny kindly 
made available to us. The three laboratory measured wave- 
lengths were also compared to the corrected wavelengths. They 
agreed to within 0.04A. There are evidently many CN lines in 
the region but none of them seems to affect the [C I] 8727 A 
line significantly. The lines in the neighbourhood of the [C I] 
line are listed in Table [j]. 



Table 1. Atomic and molecular line data for lines in the im- 
mediate neighbourhood of the 8727 A [C I] line. Only lines 
with an equivalent width W\ > 0.05mA are given. The first 
five columns give, respectively, the rest wavelength, the lower 
excitation energy, the logarithm of the product of the oscilla- 
tor strength and the statistical weight, a multiplicative factor to 
the classical van der Waals damping constant and the radiation 
damping constant. * indicates a change in A, log gf, Tq and/or 
r ra d relative to the original value in VALD 
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We also searched for telluric emission and absorption lines 
in the region of interest. These lines can influence the effective 
equivalent width of the carbon line if not properly removed. We 
found no evidence for such undesired features. 

A solar MARCS model was used to produce a synthetic 
spectrum of the Sun. This spectrum was compared with a solar 
high-resolution flux spectrum (Kurucz et al. 1984) in order to 
calibrate the oscillator strengths of some larger and more im- 
portant lines including our [C I] line (see Table [j]). The differ- 
ential analysis is consequently performed relative to the Sun. 

The carbon line is supposed to be blended by a weak Fe I 
line at 8727. 13A. In our analysis, the logy / value (log gfp c : = 
—3.93) from VALD was used. An upper limit of log gfp c T = 
—3.60 can be estimated from oscillator strengths of Fe I lines of 
the same multiplet observed in the solar spectrum (Lambert & 
Ries 1977). Using this upper limit, we find that the contribution 
from the Fe I line to the total equivalent width is less than 15% 
for the coldest and most iron rich of our programme stars, and 
less than that for the others. Also, due to the differential analy- 
sis, the uncertainty in the log gfp c T does not affect the derived 
carbon abundances by more than 10% for the hottest stars in 
the sample. 

In order to verify the assumption that the carbon line is 
responsible for most of the equivalent width of the line at 
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Table 2. Column 2-4 give the atmospheric parameters for the stars adopted from EAGLNT (Edvardsson et al. 1993) and column 
5 contains our measured carbon abundances. Column 6 shows the oxygen abundance taken from EAGLNT (they do not give an 
oxygen abundance for all the stars in our sample) and column 7 is the logarithm of the age also adopted from EAGLNT. Column 
8 and 9 contain data from the HIPPARCOS catalogue, the parallax in mas resp. the V-magnitude. Column 10 give the estimate 
of the mass from evolutionary tracks as described in the text and column 1 1 give the difference between the photometric log g 
adopted from EAGLNT and log g calculated using HIPPARCOS data 
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8727A, we studied the thermal broadening by measuring the 
FWHM of the feature in the solar flux spectrum (Kurucz et 
al. 1984). We determined the width of the combination of the 
instrumental profile and macroturbulence by fitting the neigh- 
bouring Si I lines and next calculated the 8727A feature, al- 
ternatively assuming it to be due to [C I] and Fe I, respec- 
tively. The FWHM of the line feature in the spectrum was 
measured to FWHM 872 7 feature = 0.20A while the carbon 
and the iron line in the calculated spectra were measured to 



FWHM [C ,] = 0.19A and FWHM Fe , = 0.15A, respectively. 
This indicates that the observed feature is primarily the [C I] 
line. 

4. Resulting abundances 

The abundance determinations were made by using synthetic 
spectra, calculated for a model atmosphere, tailored for each 
star. The synthetic spectra were convolved with macroturbu- 
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Table 2. continued. * indicates a derived upper limit in [C/H] 



5 



Star 




logs' 


[Fe/H] 


[C/H] 


[O/H] 


logr 9 


7T 


V 


M/M Q 


A logg 


Id. 


[K] 


[cgs] 








[Gyr] 


mas 


mag 




[cgs] 


HR 4989 


6314 


4.25 


-0.28 


-0.10 


-0.25 


0.64 


55.49 


4.90 


1.12 


0.07 


HR5338 


6177 


3.94 


-0.11 


-0.03 


0.00 


0.48 


46.74 


4.07 


1.33 


0.20 


HR 5459 


5720 


4.27 


0.15 


0.22 


0.15 


0.62 


742.12 


-0.01 


1.00 


0.04 


HR 5542 


6001 


4.09 


0.13 


0.15 


0.03 


0.69 


24.99 


6.30 


1.33 


0.06 


HR 5698 


6341 


4.04 


0.01 


0.14 


-0.02 


0.51 


29.27 


4.99 


1.58 


0.22 


HR 5723 


6532 


3.93 


-0.13 


-0.02 




0.41 


30.90 


4.92 


1.41 


0.09 


HR 5996 


5831 


4.03 


0.23 


0.17 


0.09 


0.72 


34.60 


6.32 


1.18 


-0.18 


HR6189 


6200 


3.98 


-0.56 


-0.48 


-0.48 


0.74 


23.02 


6.33 


1.18 


0.00 


HR 6243 


6435 


3.80 


0.00 


0.12 




0.28 


27.04 


4.64 


1.65 


0.14 


HR 6409 


6480 


3.94 


0.09 


0.14 


-0.04 


0.38 


19.80 


5.53 


1.58 


0.20 


HR 6569 


6675 


4.15 


-0.27 


-0.06 




0.33 


45.72 


4.76 


1.33 


0.02 


HR 6649 


6088 


4.24 


-0.34 


-0.21 


-0.29 


0.86 


30.55 


6.19 


1.05 


0.16 


HR 6907 


6382 


4.15 


0.13 


0.19 


0.01 


0.35 


27.53 


5.90 


1.41 


0.06 


HR7126 


6655 


4.16 


0.21 


0.10 




0.19 


31.24 


5.56 


1.50 


-0.01 


HR 7560 


6146 


4.14 


0.09 


0.02 




0.63 


51.57 


5.12 


1.25 


-0.07 


HR 7766 


5964 


4.22 


-0.36 


-0.19 


-0.26 


0.94 


30.84 


6.26 


1.00 


0.16 


HR 7875 


5991 


4.09 


-0.44 


-0.30 


-0.30 


0.87 


41.33 


5.11 


1.12 


0.18 


HR 8027 


6285 


4.20 


-0.37 


-0.35 


-0.31 


0.73 


35.07 


5.76 


1.12 


0.08 


HR 8041 


5813 


4.20 


0.11 


0.27 




0.89 


37.80 


6.21 


1.05 


0.01 


HR 8077 


6166 


4.05 


-0.07 


0.07 




0.62 


27.06 


5.94 


1.18 


0.10 


HR8181 


6139 


4.34 


-0.67 


-0.62 


-0.57 


0.96 


108.50 


4.21 


0.80 


0.05 


HR 8665 


6228 


4.15 


-0.32 


-0.17 




0.72 


61.54 


4.20 


1.18 


0.16 


HR 8697 


6288 


3.97 


-0.25 


-0.18 




0.56 


37.25 


5.16 


1.18 


0.01 


HR 8969 


6255 


4.16 


-0.17 


-0.10 




0.64 


72.51 


4.13 


1.12 


0.07 


HD 6434 


5813 


4.42 


-0.54 


-0.33 


-0.34 


1.10 


24.80 


7.72 


0.90 


0.06 


HD 14938 


6164 


4.09 


-0.37 


-0.30 


-0.17 


0.67 


18.58 


7.20 


1.05 


0.01 


HD 17548 


5977 


4.27 


-0.59 


-0.35* 


-0.48 


1.06 


18.90 


8.16 


1.00 


-0.13 


HD 25704 


5844 


4.43 


-0.85 


-0.68 


-0.57 


1.22 


19.02 


8.11 


0.70 


0.25 


HD 51929 


5845 


4.28 


-0.64 


-0.58 


-0.43 


1.14 


26.58 


7.39 


0.70 


0.09 


HD 69611 


5795 


4.29 


-0.58 


-0.40 


-0.28 


1.08 


20.46 


7.74 


0.90 


0.10 


HD 78747 


5824 


4.45 


-0.64 


-0.56 


-0.38 


1.08 


25.16 


7.72 


0.70 


0.19 


HD 130551 


6237 


4.25 


-0.62 


-0.36 


-0.46 


0.85 


20.94 


7.16 


0.80 


0.18 


HD 165401 


5758 


4.31 


-0.47 


-0.34 




1.18 


41.00 


6.80 


0.90 


-0.10 


HD 188815 


6181 


4.29 


-0.58 


-0.40 


-0.42 


0.86 


17.82 


7.47 


1.05 


0.13 


HD 199289 


5894 


4.38 


-1.03 


-0.66 


-0.72 


1.26 


18.94 


8.28 


0.70 


0.11 


HD 201891 


5867 


4.46 


-1.06 


-0.55* 


-0.33 


1.23 


28.26 


7.37 


0.70 


0.22 


HD 205294 


6236 


4.01 


-0.36 


-0.24 


-0.28 


0.64 


18.78 


6.86 


1.18 


-0.02 


HD 210752 


5910 


4.25 


-0.64 


-0.50* 


-0.56 


1.15 


26.57 


7.44 


0.70 


0.02 


HD 215257 


5983 


4.37 


-0.65 


-0.55 


-0.52 


1.06 


23.66 


7.41 


0.70 


0.23 


HD 218504 


5945 


4.20 


-0.62 


-0.40 


-0.43 


1.03 


16.20 


8.11 


0.70 


0.12 



lence, rotational and instrumental broadening profiles, in order 
to match the observed spectral lines. In most cases the rota- 
tional convolution was not necessary. We based the convolution 
parameters for each spectrum on the profiles of the Si I line at 
8728. OlA and the Si I line at 8742.45A. There is, however, 
no unique way to convolve the spectra to obtain the observed 
line profiles. By changing the parameters of the convolution 
in different ways we could see that the derived carbon abun- 
dances were not sensitive to the choice of convolution profile. 
The carbon abundance was next changed until a good fit to the 
observed spectrum was provided. 

The derived abundances are presented in Table |[ In the 
spectra for the stars HD 17548, HD 201891 and HD 210752 



the carbon line was too weak to allow an accurate measurement 
of the carbon abundance, so we could only determine an upper 
limit. Those stars are excluded from further analysis. 

5. Error analysis 

5.7. Errors in observed spectra and continuum location 

The largest statistical errors are introduced when measuring the 
abundance - that is when determining the continuum and fitting 
the synthetic spectrum to the observed one. These errors vary 
from star to star depending on the quality of the spectra. The 
estimated mean error for the more metal poor half of the stars 
is ±0.07 dex and for the metal rich half ±0.04 dex. 
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Table 3. Effects on logarithmic abundances derived when changing the fundamental parameters of the model at- 
mospheres. Three stars of different parameters are presented. The parameters (T e s, logg, [Fe/H], £ t ) are, respectively: 
(5965, 4.37, -0.82, 1.39), (6480, 3.94, 0.09, 2.36) and (6517, 4.01, -0.24, 2.15) 





HR 3578 




HR 6409 




HR368 




Change 


A[C/H] 


A[C/Fc] 


A[C/H] 


A[C/Fe] 


A[C/H] 


A[C/Fc] 


AT cS = +100 K 


+0.03 


-0.03 


+0.03 


-0.03 


+0.04 


-0.02 


Alogff = +0.2 


+0.07 


+0.08 


+0.07 


+0.08 


+0.06 


+0.07 


A [Fe/H] = +0.1 


+0.02 


+0.02 


+0.01 


+0.01 


+0.01 


+0.01 



5.2. Errors in fundamental parameters 

The stars in our sample are chosen among the stars analysed in 
EAGLNT, where an extensive error analysis is also provided. 
A discussion on the model atmospheres, used also in this study, 
and the derivation of the fundamental atmospheric parameters 
used for the model atmospheres are presented in EAGLNT. 

The statistical errors of the effective temperature are esti- 
mated in EAGLNT to be of the order of 50 K. In addition to 
these random errors, systematic shifts in the temperature scale 
of 50 — 100 K should be considered. This limit on systematic 
errors is consistent with the tests performed for some of the 
stars from EAGLNT by Tomkin et al. ( |1995[ ). 

The logarithm of the gravity has estimated statistical errors 
of 0.07 dex, but with possible systematic errors of 0.2 dex or 
even more, cf. the discussion below. 

The statistical errors in the model metallicity and in the mi- 
croturbulence, £ t , are estimated to be unimportant for the de- 
rived carbon and iron abundances, cf. Andersson & Edvardsson 



fl!994| ). 

The propagation of these parameter errors to the carbon 
abundance estimate has been investigated by changing the fun- 
damental parameters with a representative amount and then 
running the synthetic spectrum program to obtain the change 
in carbon abundance. This is shown in Table || for three stars, 
spanning the parameter space. The statistical error in the car- 
bon to hydrogen ratio due to the uncertainties in the funda- 
mental parameters is around 0.03 dex. The systematic error is 
larger, 0.08 dex, mainly reflecting the significance of the errors 
in logg. 

5.3. Errors due to departures from LTE and oversimplified 
models 

The carbon line being used here (2p x S — 2p x D) is forbidden 
due to the change of the orbital angular momentum quantum 
number by two units for transitions between the two states. 
These low excitation states are the first and second excited 



states of carbon. Stiirenburg & Holweger (1990) found for the 



Sun and Vega that non-LTE abundance effects of neutral carbon 
increase with the strength of a spectral line. For the Sun, the 
departure coefficients of the levels included in the [C I] 8727A 
line are negligible for all optical depths. The lowest three lev- 
els of C I are strongly coupled to the ground state through col- 
lisions. In their model for Vega (AO V) these low lying states 



show strong under-population due to the stronger UV fie ld, bu t 
are still strongly coupled to each other. Rentzsch-Holm ( 1996J ) 
showed in a non-LTE analysis of carbon that the non-LTE ef- 
fects depend on equivalent widths of the lines. In her 7000 K 
model the lowest three levels are in LTE all the way out in the 
photosphere. Non-LTE effects increase with effective tempera- 
ture. It seems therefore unlikely that the [C I] 8727A line will 
be subject to non-LTE effects. Thus, no serious systematic er- 
rors are to be expected due to our assumption of LTE. 

Other systematic errors in the model atmospheres, e.g. in 
the temperature structure, errors in the basic assumptions of 
the modelling of the atmospheres such as that of plane-parallel 
stratification and effects of unresolved binary stars, are in gen- 
eral terms discussed in EAGLNT. For the carbon abundances 
derived here, these errors may be judged to be less significant 
than abundance determinations for most elements, in view of 
the relatively low excitation energy of the line, its weakness 
and the fact that it represents the dominating ionization stage 
of carbon. The error in the carbon abundance due to the errors 
in the determination of the oscillator strength of the carbon line 
is only ±0.01 dex (Andersson & Edvardsson 1994). 



6. Relations between abundances and comparison with 
other results 

In Fig. | we show [C/H] vs. [Fe/H] for the stars observed. Ob- 
viously there is a tight correlation. Assuming a linear model, 
the spread in the points is found to be consistent with errors 
in [Fc/H] ~ 0.05 dex, as found by EAGLNT, and with errors 
in [C/H] according to the random errors estimated in the dis- 
cussion above. The slope A[C/H]/A[Fe/H] is (0.85 ± 0.03), 
with the error representing the statistical errors. What are the 
effects of the systematic errors discussed above on this rela- 
tion? A systematic trend is found when comparing the log g- 
values derived from Stromgren photometry (Edvardsson et al. 
1993) and the values obtained using parallaxes measured by the 
Hipparcos satellite. For stars of low metallicities an even lower 
abundance is derived - that is, a steeper increase of [C/H] vs. 
[Fe/H] is found - if Hipparcos parallaxes are used. The slope 
is in the latter case A [C/H] /A [Fe/H] = (0.88 ± 0.03). 

In Fig. |]a, we have plotted the [C/Fe] ratios for the ob- 
served stars as a function of [Fe/H]. There is a well defined 
relation with a small scatter around it: 

[C/Fe] = (-0.17 ± 0.03) x [Fe/H] + (0.065 ± 0.008). (3) 



B. Gustafsson et al.: The origin of carbon 



7 



0.4- 
0.2- 


-0.2- 
-0.4- 
-0.6- 
-0.8- 



***** 4 

*5* 



-0.6 -0.4 -0.2 
[Fe/H] 



0.2 



Fig. 3. [C/H] as a function of [Fe/H]. The triangles denote up- 
per limits in [C/H] and will be excluded in the analysis. A typ- 
ical error bar is shown in the lower right corner 

The standard deviation is 0.06 which again is consistent with 
the estimated random errors in carbon and iron abundances. 
Also in this case the effects on the slope from possible system- 
atic errors in log g are minute. 

The small scatter around the linear relations discussed here 
indicates that the estimates made above of random errors in the 
carbon abundances are not underestimated. The slope in the 
second relation, —0.17 ± 0.03 is significantly less steep than 
that in the relation between [O/Fe] and [Fe/H], 

[O/Fc] = (-0.36 ± 0.02) x [Fc/H] - (0.044 ± 0.010), (4) 

as obtained from the results of EAGLNT and also shown in Fig. 
Qb. The slopes as such reflect the slowly growing significance 
of iron-producing supernovae Type la during the evolution of 
the Galaxy, as compared with the rapid production of oxygen 
by Supernovae Type II (see, e.g., EAGLNT). However, the dif- 
ference between these two slopes will be the major issue in our 
discussion. 

Carbon abundances of Galactic Disk stars have been deter- 
mined by several other groups. We have 20 stars in common 
with the study of Laird ( 1985[ ) and 29 in common with Tomkin 
et al. (1995 ). Furthermore, we have a few stars in common with 



the studies of Friel & Boesgaard (1990, 1992) and Tomkin & 



Lambert (1984 ) and Clegg et al. (198 1 ). A comparison between 



the different results is shown in Fig. 5 



Andersson & Edvardsson (1994) found a slope in C/Fe, 
based on the same line as we have used, however, for fewer 
stars and lower S/N, of —0.2. They also extracted data from 
previous analyses and separated abundances based on sepa- 
rate features for comparison. The abundances found by Laird 
( |l 985 ) from CH molecular lines show a slope in the [C/Fe] vs. 
[Fc/H] of -0.10 if only stars with [Fe/H] > -1.0 are con- 
sidered. The data show considerable scatter, however. Clegg et 
al. (1981) found a mean slope of ss —0.16 ± 0.08 from highly 
excited C I lines, the [C I] 8727A line and CH lines in F and G 
main-sequence stars. Their data demonstrate systematic differ- 
ences between results from the different lines used: The highly 
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Fig. 4. a [C/Fe] as a function of [Fe/H]. The asterisks denote 
abundances of observed stars and the full line is the determined 
slope given in Eq. 3. b [O/Fe] as a function of [Fe/H] with data 
taken from Edvardsson et al. (1993). The slope calculated from 
this data is given in Eq. 4. Typical error bars are shown in the 
lower left corners 



excited C I and the [C I] 8727A line give a slope of -0.18 and 
the CH-based abundances give —0.27, while a slope of 0.00 
was obtained from C2 lines. Friel & Boesgaard (1990, 1992) 
show from highly excited C I lines in fi eld sta rs a slope of 
—0.18, while , more recently, Tomkin et al. ( 1995 ) derived from 
the C I line at 7100A of field disk stars a slope of -0.28. 

In Fig. H we also note a small offset of our data; e.g. relative 
to those of Tomkin et al. (1995) we tend to obtain [C/H] ratios 
greater by 0.1 dex. This may be due to the use of different car- 
bon lines. If our results were true they also indicate a tendency 
for the Sun to be somewhat carbon poor or iron rich (cf. also 
Fig. ||). This latter circumstance should be compared with the 
possible tendency for the Sun to be somewhat poor also in other 
elements like Mg, Al and Si relative to Fe (Gustafsson 1998). 
In view of the possible systematic errors involved, e.g., in the 
gravity scale, we do not consider our offsets from the identity 
line in Fig. || and our finding for the solar type stars to have 
non-solar [C/Fe] ratios to be very significant. This, however, 
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Fig. 5. Comparison of derived [C /H] for the stars in common 
between earlier studies and present study. Plusses (+) are from 
Clegg et al. (1981), crosses (x) are from Tomkin & Lam- 
bert (1984), points (•) are from Laird (1985), squares (□) are 
from Friel & Boesgaard (1990), diamonds (O) are from Friel 
& Boesgaard (1992), and open circles (o) are from Tomkin et 
al. (1995). The dashed line shows the loci of identical results 

deserves further study. Our data show a smaller spread in the 
[C/Fe] vs. [Fe/H] relation than do the other studies. We ascribe 
this to the high quality of our data, and the insensitivity of the 
8727A [C I] line to temperature. 

7. Discussion 

7.1. Possible sources of carbon 

As was mentioned in the Introduction a number of different 
production sites have been suggested for carbon: supernovae, 
novae, Wolf Rayet stars, intermediate and low mass stars in 
connection with the planetary-nebulae stages or even before by 
superwinds at the end of the red-giant phase. Here, we shall 
discuss these different possibilities in the light of our observa- 
tions. 

If carbon were produced essentially by SNe II, with yields 
relatively independent of the initial metal abundances of the 
massive stars from which the supernovae originate as is the 
case for oxygen, then one would expect the slope in the two 
relations in Fig. 4 to be the same. This is because the yield 
of carbon produced by a generation of massive stars would 
scale with that of oxygen. True enough, the C/O ratios of 
the nuclei produced and expelled by SNe II of different mass 
may be different. E.g., from the tables of Nomoto et al. (1997) 
we find 12 C/ 16 ratios varying from 0.23 (at 15M ), 0.21 
(18.M ), 0.049 (25M Q ) and 0.016 (407W Q ). An analogous 
trend is found from the grid of supernovae models of Woosley 
& Weaver ( |1995j ) (henceforth referred to as WW95) although 
their C/O ratios are greater, reflecting their lower rates for the 
critical 12 C(a, 7) 16 reaction. However, the time scale differ- 
ence between stars in the lower end of the mass interval that 
produce SNe II and stars of very great masses is still much 
smaller than the characteristic time scale for SNe la, produc- 



ing iron. Therefore, these different yields of supernovae of dif- 
ferent mass are not expected to lead to any significant slope 
differences in Fig. 4. If the C/O ratios were dependent on the 
initial metallicity of the star evolving to a supernova this would, 
however, lead to a difference. In effect, if the most metal-poor 
stars would produce relatively more oxygen one could at least 
qualitatively reproduce the slopes in Fig. ||. However, the SN 
II models with different metallicity do not predict any substan- 
tial differences in the C/O ratios of the yields. E.g., the mod- 
els with solar initial abundances of WW95 have yields with 
O/C = 7.1, when integrated with a mass distribution corre- 
sponding to a Salpeter IMF, while those for 1/100 solar metal- 
licity lead to O/C = 6.6. Also note that the absolute amount 
of carbon produced in the models mentioned, as compared with 
other SN II products, e.g., Mg, is not sufficient to explain the 
abundance difference of a factor of 10 found, e.g. in the Sun, 



cf. WW95 and Thielemann et al. ( 1996 ). Although there are 
uncertainties still in predicted carbon yields for SNe II due to 
the uncertainty in the 12 C(a, 7) 16 rate and the details in the 
treatment of convection, cf. WW95, the uncertainties do not 
seem large enough to admit an increase of the carbon yields by 
a factor of 5. Obviously, another source of carbon is needed. 

Supernovae of Type la, believed to produce most of the 
iron, do not contribute significant amounts of carbon. E.g., 



the various models listed by Nomoto et al. (1997) give carbon 
yields smaller than 10% of their iron yields. 

Novae may also be ruled out as important sites for the pro- 
duction of carbon since the frequency of such events and their 
mass loss are too low to enrich the interstellar medium enough 
with common elements such as 12 C (Gehrz et al. 1998). 

Massive stars may, however, contribute processed material 
in earlier phases, due to their massive radiation-driven winds. 
These winds may then contain He and C due to hydrogen and 
helium burning, but on ly smaller amounts of O and heavier el- 
ements. Maeder ( 1992 ) suggested yields from stars more mas- 
sive than 25.M0 to be strongly dependent on metallicity, re- 
flecting the fact that more metal rich gas has a much greater 
cross section for radiation in the ultraviolet. E.g., a 40.M© solar 
composition model is suggested to eject about 10% of its mass 
as He and 10% as C into the interstellar medium, while a corre- 
sponding model with 1/20 of the solar metallicity hardly looses 
anything. These calculations have been detailed by Portinari et 
al. (1998) who also produce results for intermediate metallici- 
ties. In their models, C/O of the winds is less than solar for ini- 
tial masses in the interval 12 - 3OA40. Not until the mass is on 
the order of 40A^o or greater does C/O become greater than 
one. Although there is a general agreement between observed 
and predicted relations for [C/Fc] vs. [Fe/H] and [O/Fc] vs. 
[Fc/H], the fit is not perfect. E.g., for carbon a bump with high 
carbon abundances is predicted around [Fc/H] = —0.5, with 
a steep slope down in [C/Fe] when proceeding towards solar 
[Fe/H]. This is not seen in our observations. On these grounds, 
however, the origin of carbon in massive stars may certainly not 
be refuted. The substantial uncertainties in calculated yields, 
both from supernovae and the massive winds in earlier stages 
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of massive stars, may well explain these deviations from obser- 
vations. 

A more empirical estimate verifies the role that massive 
stars may play through their winds. The number of WR stars 
within 3kpc from the Sun, a volume for which the sample is 
thought to be complete, is about 60 (Conti 1988). Of these, 
about half are WC stars which then implies a number of WC 
stars in the Galaxy on the order of 300. Typical mass-loss rates 
from WC stars are 7 • 10 -5 A^©/year and C/He mass ratios 
range in the interval 0.5 - 3.3 (de Freitas & Machado 1988). 
This implies a total annual carbon contribution to the ISM from 
the WC stars of about 0.01.M©, which is a very considerable 
number. In model calculations to be described below, we find a 
value fairly consistent with that, c onside ring the uncertainties: 
with the yields of Portinari et al. ( 1998 ) we obtain a total car- 
bon contribution from massive stellar winds of 0.0035.Mo per 
year, for a star formation rate of lyVf© per year in the Galaxy, 
and a Salpeter IMF. 

Stars of intermediate and low masses (< lOvVl©) may 
also be significant contributors of carbon on the galactic scale, 
as was suggested by Tin sley (|1978 ), Sarmiento & Peimbert 
( 1985 ) and Timmes et al. (1995). The basic argument for these 
authors was that the estimated contribution from supernovae 
and novae was not enough to account for the observed carbon 
abundance. The theoretical yields of carbon from intermedi- 
ate and low mass stars are, however, highly uncertain, not the 
least for the numerous low-mass stars. The pioneering calcu- 
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Fig. 6. [C /H] as a function of log rg. The scatter in this relation 
is much greater than in the [C/H] vs. [Fe/H] relation (Fig. 3), 
partly due to errors in the age estimates 



the halo phase of about 1-2 Gyrs admits the contribution of 
carbon by intermediate mass stars to the initial abundance of 
the Galacti c Disk, as predicted from the yields of Renzini & 
Voli ( 198 1 ), but that the lack of a peak in the C /O ratio around 
[Fe/H] = —1.0 suggests other production sites. They conclude 
from their model calculations that an origin in massive winds 
as suggested by Maeder (1992) is most probable. 



lation of yields by Renzini & Voli (|1981|) were not success- ? 2 Tes(s ofpossMe carbon sources 



ful in predicting observed C/O and N/O ratios for stars and 
planetary nebulae, nor could they reproduce the carbon stars of 
relatively low luminosities in the Magellanic clouds (cf. Note 
added in proofs by these authors). These calculations are now 
surpassed by those of Marigo et al. (1996, 1998) who have 
published yields based on semi-analytical models of AGB stars 
with mass loss, with several free parameters to describe dredge- 
up of processed material. They find significant contributions of 
carbon from stars in the mass interval 1.5 - AM.Q. Forestini & 



Charbonnel ( 1997 ) have calculated evolutionary models to the 
tip of the AGB with pulsations included and found significant 
contributions of carbon for the lowest masses, amounting to a 
few percent of a solar mass for stars around 2 - 2.5.M©. Both 
these papers suggest greater yields of carbon for low metal- 
licities (Fig. ^|). Empirically, the finding that 50% of the plan- 
etary nebulae are carbon rich (with C/O > 1), (Zuckerman 
& Aller 1986; Rola & Stasinska 1994) and a PN birth rate of 
1 per year in the Galaxy (Pottasch 1992) and a PN mass of 
O.3A^0, suggests a carbon contribution of about 0.003A^q per 
year to the ISM. With such a contribution, the PNe could be a 
major source of carbon. One may ask whether mass loss from 
carbon stars, i.e. before the PN stage, also could be a major 
source. Adopting typical mass-loss rates (see Olofsson et al. 
1993, 1996; Jura & Kleinmann 1989), C/O ratios (Lambert 
et al. 1986) and carbon-star frequencies we find a total carbon 
contribution of 2 ■ 10 -4 A^e/year, which should not be very 
significant in the global c arbon budget of the Galaxy. We note 
also that Prantzos et al. ( 1994 ), argue that the time scale for 



One may now ask which of the two major remaining carbon 
sources - massive stars through massive stellar winds, presum- 
ably in the Wolf Rayet WC stage, or the contribution from in- 
termediate/low mass stars through planetary nebulae - is the 
most significant source. These two different sources could re- 
veal themselves in different strengths of the correlations be- 
tween carbon abundance and metallicity, and age, respectively. 
The WC origin of carbon should be expected to lead to a 
stronger correlation of the carbon abundance with metallicity, 
while if the age is the dominating factor, determining the car- 
bon abundance as in the planetary nebulae case, one would ex- 
pect a stronger correlation between [C/H] and age. 

In contrast to the strong correlation between the carbon 
abundance and metallicity, as illustrated in Fig. Bl a much 
greater scatter is found in the [C/H] vs. stellar age relation 
(Fig. ||). A similarly great scatter is found in the [Fe/H]-age 
diagramme (see EAGLNT, their Fig. 14a). This could indicate 
a cosmic spread in [Fe/H] at a given age, which then indicates 
that [C/H] is more tightly related to metallicity than to age. 
Another possibility is that the scatter in Fig. || is due to random 
errors in age. If so, these errors must be as large as 0.2 dex or 
around 60%. The estimated errors in the ages are more likely 
on the order of 25%. Therefore, it seems likely that there is a 
real spread in the metallicity-age relation and that the primary 
dependence is on metallicity. This argument is weakened, how- 
ever, by the possibility that the spread in the [Fe/H] — age rela- 
tion is due to the diffusion of stars formed at different distances 
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Fig. 7. [C/O] as a function of [Fe/H]. We can see a grad- 
ual increase in [C/O] with [Fe/H]. This observed enhance- 
ment of carbon to the ISM could be a result of carbon pro- 
duction in massive stars with metallicity dependent yields or 
intermediate- and low-mass stars, or possibly a mixture 



from the Galactic centre into orbits close to the solar one (Wie- 
len et al. 1996). 

We here suggest another test, built on the following ar- 
gument: If carbon is produced in massive stars, but with 
metal-abundance dependent yields due to the radiatively driven 
winds, one would expect a gradual increase of the C/O ratio 
with [Fe/H]. This is also observed in the Disk, (Fig. |7|). Simi- 
larly, if less massive stars produce carbon and expel it as PNe, 
one again expects a gradual increase of C/O in the Disk, re- 
flecting the difference in time scale between the carbon pro- 
ducing stars and the massive stars producing oxygen. 

Our test is built on the fact that these two factors behind 
the C/O increase - metallicity and time-scale - are interre- 
lated in different ways in different galaxies. Thus, irregular 
galaxies, like the Magellanic Clouds, have a lower metallicity 
([Fe/H]i rr ), while the time scale for a more or less continuous 
star formation history may be of the same order of magnitude 
as for the Galactic Disk. For these irregulars we therefore ex- 
pect the massive star contribution of carbon relative to oxygen 
to be be comparable to that from massive stars in the Disk, 
when its metal-abundance was [Fe/H]i rr . Thus the same slope 
should be delineated in the [C/O]- [Fe/H] diagram by the ir- 
regulars and the Galactic Disk stars. For the contribution of 
carbon from intermediate- and low-mass stars, however, we ex- 
pect a smaller contribution in the Disk when the metal abun- 
dance was [Fe/H]i rr if the latter is small enough (< —0.5), 
since time then had still not admitted these stars to contribute 
much. In the more slowly evolving irregulars, however, these 
contributions are fully developed. Thus, we predict a greater 
slope in the [C/0]-[Fe/H] diagramme for the Galactic Disk 
stars if intermediate- and low-mass stars are more significant 
as carbon sources. 

It should be noted that Edmunds & Pagel ( 1978| ) pointed 
out that, provided that nitrogen would be formed as a primary 
element in low-mass stars, the N/O-O/H diagramme could be 



used for dating galaxy populations. Here, we use a similar idea, 
for another purpose. 

7.3. Models of carbon synthesis in galaxies 

We shall now illustrate our discussion by using models of 
Galactic chemical evolution. Such models must, however, be 
used with care - in particular when compared quantitatively to 
real observations - since they are marred by the uncertainties in 
the in-going data, in particular in yields. Also, the uncertainties 
in the basic assumptions made in the models are significant. 
This is illustrated in Fig. 3 of Carigi (1994) where a number of 
different models from the literature are intercompared. 

Our chemical evolution model is based on the analytical 
model of the local outer Disk, developed by Pagel & Taut- 



vaisiene (1995). This model uses the inflow formalism of Clay- 



ton (1985) and the delayed production approximation intro- 



duced by Pagel (1989) 



The differential equations describing the chemical evolu- 
tion are: 



dgi 
du 

dZ 2 
du 



F 
tog 



Zi = y x {Z{u)) 



^Z 2 
tog 



= 



dZ 2 
du 



^ + —Z 2 = y 2 (Z(u)) 



i£u < wA; 
g(u — wA) 

g( u ) 



ifu > luA. 



(5) 
(6) 
(7) 



where the total mass fraction of a species, Z\ = Z\ + Z 2 , is 
the sum of the mass fractions from instantaneous and delayed 
production. Here, Z is the overall metallicity (per mass). The 
time-like integration variable u, is defined as 



w(t')dt', 



(8) 



where u> (= 0.3 Gyr™ 1 ) is the transition probability for gas to 
change into stars in unit time at time t. The function describing 
the inflow is of the form 



F k 
ujg u + uq 



(9) 



and can be chosen to be switched on at a specific time u\. F is 
the mass inflow rate, while the parameters k and uo are choosen 
to fix the ratio of final to initial mass of the system. In this 
model, u \ — .14, k = 3 and uq = 1.3, as in Pagel & Taut- 
vaisiene ( 1995 ). Note that before the inflow begins, the second 
term on the left hand side of Eqs. 5, 6 and 7 is identically zero. 
Eventually, the evolution of the gas mass, g(u), is given by 



dg 

du. 



F 



(10) 



All parameter values describing the chemical evolution are 
the same as those chosen by Pagel & Tautvaisiene except for 
their choice of constant net yields. We explore the effects of 
metallicity dependent yields P\m{Z{u)), where Z(u) is the 
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Fig. 8. Synthetic [C/O] ratio of the local disk and dwarf irreg- 
ular galaxies derived from a chemical evolution model using 
metallicity dependent yields (Portinari et al. 1998). The full 
line shows the [C/O] ratio of the Galactic Disk with carbon 
produced by high mass stars, the dashed line shows the Disk 
with products from intermediate and low mass stars, the dotted 
line is the combined [C/O] ratio with carbon from high, inter- 
mediate and low mass stars in the Disk. The asterisks (*) repre- 
sent dwarf irregular galaxy models with carbon from high mass 
stars, the plusses (+) are dwarf irregular galaxies with carbon 
from intermediate and low mass stars and the points (•) rep- 
resent dwarf irregular galaxies with carbon produced by high, 
intermediate and low mass stars. The sharp upturn in [C/O] 
at [Z] = —0.6 is a result of the highly increasing theoretical 
carbon yield of high mass stars with increasing metallicity 



overall metallicity. The net yields (yi in Eqs. 5 and 7) are mass- 
integrated stellar yields divided by the lock-up fraction, which 
is treated in the instantaneous recycling approximation. 

We have used carbon and oxygen yields of Portinari et al. 



(1998) for high mass stars and carbon yields of Marigo et al. 
(1996, 1998) and Forestini & Charbonnel (|1997|) for intermedi- 



ate and low mass stars. With an adopted linear time-metallicity 
relation normalized to the solar value, i.e. Z = 0.0055u, we 
can calculate the evolution of the [C/O] ratio with Z (Fig. ^|). 
Oxygen was assumed to be recycled instantaneously while car- 
bon was divided into an instantaneous part (representing car- 
bon production by high-mass stars) and several delayed pro- 
duction parts with different time delays wA, from ui A = 0.015 
to a; A = 1.000, corresponding to carbon production by stars 
with masses from 6 - 9M & to 1 - 2M Q . 

The dwarf irregulars were simulated like local Disk sys- 
tems (i.e. same parameters as the Disk), however, with differ- 
ent slopes of the time-metallicity relation in order to mimic 
the slower metal production of heavy elements leading to a 
lower present metallicity. All systems were evolved for 15 Gyrs 
{u = 4.5) and thereafter the [C/O] ratio was read off for com- 
parison with observations. 

The results of the simulations are shown in Fig. ^|. It is 
seen that the C/O ratio for the Galactic Disk, if the car- 
bon is produced by massive stars, shows a gradual and rather 
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Fig. 9. a Theoretical stellar carbon yields (in A4q) for interme- 
diate and low mass stars. The full line is the yield at Z = 0.02 
and the dashed line is the yield at Z = 0.008 from Marigo 
et al. (1996, 1998). The dash-dotted and the dotted lines show 
the yield at Z = 0.02 and Z = 0.005 respectively, calculated 
by Forestini & Charbonnel (1997). Alternative increased car- 
bon yields for low mass stars (0.95 - 2.5.Mq) are indicated as 
dotted horizontal lines, b Synthetic [C/O] ratios in the Galac- 
tic Disk with the increased carbon yields for low mass stars. 
Symbols as in Fig. 8 



steep increase with Disk star metallicities, reflecting the metal- 
abundance sensitivity of the radiatively driven winds. The mod- 
els representing dwarf irregulars line up nicely along the Disk 
curve. When carbon production from intermediate and low 
mass stars is added, the dwarf irregulars are elevated above the 
C/O ratios of the Disk, in particular for the low metallicity sys- 
tems. Thus, in this case the dwarf irregulars should delineate a 
relation with smaller slope in the diagramme, as was argued 
on qualitative grounds above. In the extreme case (lower part 
of the figure) when only intermediate- and low-mass stars con- 
tribute the slope difference is very pronounced. 

Since the relative contribution to the production of carbon 
of stars of different masses is highly uncertain we have ad hoc 
increased the yields of Marigo et al. for the stars with masses 
below 2.5A^0 according to Fig. ||a. Such a change will guar- 
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antee that the frequency of planetary nebulae with C/O > 1 
gets as high (or even higer) than observed, and that the absolute 
amount of carbon produced is as high as observed. This is il- 
lustrated in Fig. ^b. Again we notice the very pronounced slope 
difference between the Disk star relation and that of dwarf ir- 
regulars. 

7.4. Results of the tests 

Before we discuss the results of the test it should be noted that 
our observations of stars in the Solar neighbourhood may not 
map the chemical evolution all the way back to the time of for- 
mation of the Disk. Short lived intermediate-mass stars could 
have enriched the ISM, provided that they are sites of high car- 
bon production. This is in clear contradiction, though, to the 
theoretical yield calculations (Fig. ||a) which indicate that the 
intermediate-mass stars in the upper mass region do not pro- 
duce carbon in large amounts. They are also less numerous in 
comparison to stars of lower mass which further limits their rel- 
ative contribution to the ISM. Moreover, if intermediate-mass 
stars were the main origin of carbon in the Disk in general, the 
carbon yield has to be metal dependent as for the massive stars, 
again in contradiction to the theoretical yields which rather de- 
crease with increasing metallicity. However, the yields are un- 
certain and observationally, we are not able to exclude them as 
carbon producers of some importance. 

In order to carry out the test empirically we first use the 
observations of massive supergiants in the Magellanic Clouds. 
The determinations for supergiants in the SMC (Hill et al. 
1997) and LMC (Hill et al. 1995) are plotted in Fig. [R]. The re- 
lation for the Disk stars was obtained by combining Eqs. 3 and 
4. From this we find that the stars in the Magellanic Clouds 
seem to take positions in the diagramme characteristic of the 
Galactic Disk stars which does not support the hypothesis of 
C being made in intermediate- or low-mass stars. However, the 
supergiant CNO abundances may be affected by dredge-up of 
CNO processed material, which makes the value of this test 
uncertain. 

In Fig. we have also plotted the location of H II re- 
gions in various metal-poor galaxies according to Kobulnicky 



& Skillman (1998) and references therein and the locus of halo 



dwarf stars according to Tomkin et al. (1992). There seems to 



be a tendency for the halo stars to have low C/O ratios as com- 
pared with Disk stars. This may, however, be due to systematic 
errors is the analysis of the halo stars, or possibly of the Disk 
stars (the latter is suggested by the deviating position of the sun 
in Fig. 10), but the suggested tendency is worth further investi- 
gation. 

In Fig. [l(] we find the galaxies to define a slope that seems 
as steep or even steeper than that of the Galactic Disk stars. Al- 
though adopting an [O/Fe] slope of —0.40 (cf. Eq. 4 and Fig. 
Qb) the picture does not alter significantly. Actually, if we drop 
the mathematically derived equations (Eqs. 3 and 4), a some- 
what higher value of both the [O/Fe]- and [C/Fe] slope may 
be argued for just by inspecting Fig. Q by eye. The C/O slope 
would then again decrease towards the adopted value (Fig. |l0| ). 
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Fig. 10. Observed C/O ratios in various metal-poor galaxies. 
Open circles (o) represent data for H II regions in different 
galaxies from Kobulnicky & Skillman (1998) and references 
therein, the square (□) is the mean log C/O ratio calculated 
from several observations of the Orion nebula. The triangles 
(A) are C/O ratios deduced from supergiant stars in LMC and 
SMC by Hill et al. (1995, 1997). The full line indicates the cal- 
culated slope of the Galactic Disk stars, derived from observed 
carbon and oxygen abundances. The hatched area in the lower 
part of the figure indicates the location of halo dwarf stars, stud- 
ied by Tomkin et al. (1992). The location in the horizontal di- 
rection of these stars is relatively uncertain due to temperature 
uncertainties. The Sun is also indicated in the figure as a 0. 
The slope of the Disk stars is not steeper than the slope of the 
dwarf irregulars, indicating that low-mass stars are not signifi- 
cant carbon contributors relative to metal-rich high-mass stars. 
The role of the short lived intermediate-mass stars should be 
further investigated 



Even if the slope of the Disk stars is further steepened when 
the halo stars are added, the locii of the metal-poor galaxies 
do not seem compatible with the hypothesis that carbon has 
been formed mainly in stars of lower mass. The most direct 
conclusion from this would thus be that the production of car- 
bon occurs at a time-scale characteristic of high-mass or possi- 
bly intermediate-mass stars of relatively high mass, but that the 
yields are metal-abundance dependent, just as is expected from 
massive stellar radiatively driven winds, e.g. from WC stars. 

There may be several arguments against a test of this char- 
acter. One may be the possibility that carbon, and oxygen as 
well as iron, could be severely depleted by grain formation 
in the interstellar medium of the irregular galaxies. A deple- 
tion independent of metallicity for carbon, but increasing with 
metallicity for oxygen, seems possible and would qualitatively 
lead to the effect observed. It is, however, questionable whether 
this effect could be as great as required to explain the full effect 
or to fully compensate for the effect due to carbon formation by 
low-mass stars. 

Another possibility may be that the remnants from the oxy- 
gen producing SNe are lost from the irregulars, and more so 
later in their evolution when the gas density in the galaxy and 
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its surroundings is less, while the slower carbon-rich winds 
from the PNe are mostly retained. This would, however, dimin- 
ish the slope of the irregulars in the diagramme, which would 
then strengthen our arguments for the high-mass stars as main 
contributors of carbon. 

The dwarf irregular galaxies could be systems of much 
smaller age than the Galactic Disk. E.g., even if Garnett et al. 
( 1997 ) have pointed out that a population of stars older than the 
present star burst is probable in I Zw 18, the mean age of stars 
in this system could still be much younger than the Galactic 
Disk. However, at least in the Large Magellanic Cloud there is 
strong evidence for a dominating population of an age compa- 
rable to that of the Galactic Disk (Ardeberg et al. 1997). 

Finally, if the main source of carbon is high-mass stars one 
has to explain why planetary nebulae, found above to provide 
as much as 0. 003.M© of carbon per year to the Galaxy today, 
do not constitute a major carbon source. The statistics of PNe 
must comprize low-mass stars (< 2A4q) as a major group; in 
view of the great PN birth rate. Therefore, if the masses, the 
birth rate or the carbon abundances of the PNe are not seri- 
ously overestimated, one has to advocate metal-abundance sen- 
sitive yields for these, such that the amount of carbon produced 
be greater the higher the metal-abundance. This is contrary to 
the tendency found in model calculations both by M arigo et al. 
(1996, 1998) and by Forestini & Charbonnel (|1997|). 



8. Conclusions 

We have determined accurate carbon abundances for a set 
of 80 Galactic Disk stars. We find the resulting [C/H] val- 
ues to correlate very well with the iron abundances [Fc/H], 
with a slope, however, that is smaller than unity such that 
A[C/Fe]/A[Fe/H] = -0.17 ± 0.03. From this, and the cor- 
responding variation of [C/O], we conclude by comparing to 
dwarf irregular galaxies of different metallicity that carbon 
is predominantly produced in the Galactic Disk by massive 
stars, ejecting carbon in radiatively driven massive winds, in 
particular during their Wolf Rayet stage. Estimates of yields 
from massive stars support this conclusion, even if yields from 
intermediate- and low-mass stars also suggest those to be of 
some significance as sources of carbon. A particular problem 
is raised by the estimates that about half the planetary nebulae 
may be carbon rich (with C/O > 1); if so one has to explain 
where this excess carbon has gone. In particular, why is it not 
visible as a steeper slope for the Disk stars than for the dwarf 
galaxies in the [C/O] - [O/H] diagramme? 

The systematic study of carbon abundances of WR stars, 
carbon stars and planetary nebulae of different progenitor mass 
and metallicity in the Galaxy and neighbouring galaxies should 
be persued, in order to establish more reliable empirical yields 
for high-, intermediate- and low-mass stars. These abundances 
should also be related to the abundances of nitrogen and the s 
elements which are thought to be formed in intermediate- and 
low-mass stars. A test similar to that applied here with metal- 
poor galaxies should be possible and rewarding by using radial 
gradients in C/O and O/H measured for disk galaxies. Also, 



isotopic characteristics of WR stars in solar-system material, 
e.g. as presolar grains in meteorites, should be further looked 
for. At present that situation is not clear; according to the recent 
review by Arnould et al. (1997) there is to date no clear and un- 
ambiguous signature of the existence in meteorites of presolar 
grains of WR origin identified in the laboratory, while grains of 
AGB star origin have most probably been identified (Anders & 
Zinner 1993). 
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